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regions.  Media  used  for  these  studies  included  high  purity  water,  tap  water, 
physiological  (0.9%)  saline  solution,  and  bovine  vitreous.  10ns  pulses  at  532  nm 
and  60  ps  and  90  fs  pulses  at  580  nm  were  focused  into  a  sample  to  produce  LIB. 
Probit  analysis  was  used  to  determine  the  50%  probability  threshold  value  (ED^q)  as 
the  temperature  of  the  media  was  varied.  Additional  data  was  obtained  by  keeping 
pulse  energy  constant  and  varying  temperature.  ^D  _  values  for  LIB  showed  no 
consistent  dependence  on  the  temperature  of  the  medium.  The  theory  of  the  tempera¬ 
ture  dependence  of  LIB  and  experimental  observations  for  all  pulse  durations  and 
their  implications  for  retinal  damage  are  discussed. 
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ABSTRACT 

The  physical  properties  of  laser-induced  optical  breakdown  (LIB)  in  highly  transparent, 
dispersive  media,  like  that  found  in  the  eye,  are  of  great  interest  to  the  ophthalmic  community.  We 
examined  the  temperature  dependent  characteristics  of  LIB  thresholds  in  media  with  a  temperature  range 
of  20°C  to  80“C  using  nanos'*<''' ,d,  picosecond,  and  femtosecond  pulses  produced  in  the  visible  and  near 
infrared  spectral  regions  \  '  i'  used  for  these  studies  included  high  purity  water,  tap  water, 

physiological  (0.9%)  salin  .uO  „  and  bovine  vitreous'.  10  ns  pulses  at  532  nm  and  60  ps  and  90  fs 
pulses  at  580  nm  were  focused  <.o  a  sample  to  produce  LIB.  Probit  analysis  was  used  to  determine  the 
50%  probability  threshold  value  (ED^  as  the  temperature  of  the  media  was  varied.  Additional  data  was 
obtained  by  keeping  pulse  energy  constant  and  varying  temperature.  EDjq  values  for  LIB  showed  no 
consistent  dependence  on  the  temperature  of  the  medium.  The  theory  of  the  temperature  dependence  of 
LIB  and  the  experimental  observations  for  all  pulsu  durations  and  their  implications  for  retinal  damage 
are  discussed. 


2.  INTRODUCTION 

The  practical  application  of  laser  induced  breakdown  (LIB)  has  been  studied  extensively  in  the 
medical  community^.  The  characterization  of  the  propagation  of  ultrashort  (subnanosecond)  laser  pulses 
through  various  media  and  the  interaction  of  the  pulse  with  the  media  is  of  fundamental  importance  in 
these  studies.  To  the  medical  community,  the  generation  and  control  of  the  plasmas  produced  by  LIB  has 
direct  application  to  ophthalmic  surgery  as  a  tool  for  ocular  tissue  disruption.  Currently,  LIB  is  used  in 
ophthalmic  surgery  in  three  different  techniques:  to  remove  the  opaque  posterior  capsule  (capsulotomy)  in 
cataract  surgery  after  the  intraocular  lens  has  been  implanted^;  to  remove  pressure  in  the  eye  by  piercing 
the  iris  in  patients  with  glaucoma^;  and  to  remove  vitreal  strands  that  form  in  the  anterior  of  the  vitreous^. 
LIB  and  its  associated  bubble  formation  and  acoustic  shock  wave  generation  may  also  cause  hemorrhagic 
lesions  when  laser  light  is  focused  on  the  retina*.  The  characterization  of  LIB,  in  particular  an  accurate 
measure  of  the  threshold,  remains  a  top  priority  in  order  to  study  the  true  interaction  of  laser  radiation 
with  the  eye. 

LIB  occurs  when  laser  irradiance  reaches  a  threshold  for  a  particular  medium.  This  threshold 
depends  on  numerous  pulse  characteristics,  such  as  duration,  wavelength,  focal  spot  size,  spatial/temporal 
profile,  and  medium  characteristics  such  as  impurity  concentration^.  Thresholds  have  been  determined  in 
the  past  for  various  focal  spot  sizes^,  and  different  focal  length  lenses’,  but  often  lack  detailed  supporting 
information  about  the  pulse  spatial/temporal  profile  and  the  medium's  physical  characteristics.  Our  goal, 
therefore,  was  to  quantitatively  measure  one  medium  characteristic  that  has  heretofore  been  ignored. 

When  a  sufficient  amount  of  thermal  energy  (kaT,  ka  -  Boltzmann's  constant)  is  introduced  into 
the  sample  not  only  will  the  population  of  translational,  rotational,  and  vibrational  excited  states  increase, 
but  the  population  of  higher  electronic  states  as  well.  We  applied  this  idea  to  LIB  phenomena  observed 
in  physiological  solutions  like  tap  water,  saline,  and  bovine  vitreous.  We  assume  that  as  the  temperature 
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increased,  the  population  of  electrons  residing  in  excited  states  would  increase.  With  an  increase  of 
excited  electrons  within  the  Rayleigh  range  of  the  focused  pulse,  statistically  the  number  of  sufficiently 
energetic  electrons  available  to  interact  with  the  pulse  (photons)  and  initiate  a  breakdown  event  will 
increase  with  the  rise  in  temperature.  Providing  our  assumptions  held  true,  we  hypothesize  a  lowering  of 
the  threshold  energy  (or  irradiance)  required  to  induce  breakdown  with  increased  sample  temperature. 

We  believe  that  if  an  accurate  measure  of  a  temperature  dependence  could  be  made  in  the  range 
of  20°C  to  80°C,  then  a  model  could  be  developed  to  determine  the  activation  energy  required  for  LIB 
initiation  and  potentially  determine  the  species  most  likely  to  contribute  the  initiating  electrons.  Further, 
and  more  importantly  to  the  study  of  ocular  damage  from  subnanosecond  laser  pulses,  most  threshold 
studies  of  LIB  associated  with  safety  or  medical  research  report  data  obtained  at  21  -  25  °C  instead  of  35  - 
37  °C.  We  seek  to  determine  if  the  difference  in  temperature  makes  any  significant  contribution  to  the 
amount  of  energy  per  subnanosecond  pulse  required  to  initiate  LIB  in  vivo. 

3.  LASER  INDUCED  BREAKDOWN 

LIB  is  an  all-or-none  process,  yet  it  is  highly  statistical  and  probabilistic  and  hence,  threshold 
measurements  tend  to  be  inexact.  Thresholds  typically  are  represented  as  an  energy  or  an  irradiance  at 
which  there  is  a  50%  or  100%  probability  that  LIB  will  occur.  Thresholds  also  have  a  large  standard 
deviation.  LIB  occurs  by  either  avalanche  or  multiphoton  ionization.  Avalanche  ionization  occurs  when 
a  seed  electron  absorbs  photons,  accelerates,  and  frees  bound  electrons  by  collisional  ionization.  These 
electrons  in  turn  accelerate  and  free  more  electrons,  and  the  process  cascades  into  a  volume  of  plasma. 
Avalanche  ionization  is  likely  to  occur  in  the  nanosecond  time  regime  where  the  pulse  is  relatively  long 
compared  to  the  collisional  ionization  dme.  This  mechanism  is  highly  impurity  dependent.  Multiphoton 
ionization  occurs  when  each  bound  electron  is  ionized  by  absorption  of  multiple  photons.  Each  electron 
acts  independently  and  no  electron-to-electron  interactions  are  required.  Multiphoton  ionization  is  likely 
to  occur  in  the  femtosecond  time  regime  and  a  transition  between  the  two  processes  occurs  in  the 
picosecond  regime.  With  multiphoton  ionization,  electrons  are  stripped  and  ions  are  created 
simultaneously  with  very  little  impurity  concentration  dependence'®  ”. 

4.  EXPERIMENTAL  SETUP 


4.1.  Laser  System  Setup 

The  ultrashort  laser  system  setup  for  producing  nanosecond  to  femtosecond  pulses  is  shown  in 
Figure  1.  The  system  is  made  from  off-the-shelf  components  built  by  Spectra-Physics.  The  primary  pulse 
comes  from  a  Model  3800  Mode-locked  Nd;YAG  laser  which  generates  80  ps  pulses  at  82  MHz  pulse 
repetition  frequency  (PRF),  1064  nm,  and  approximately  1 1 .5  W.  The  pulses  are  compressed  to  5  ps  with 
a  fiber/grating  compressor  and  frequency  doubled  to  532  nm  in  the  Model  3695  Optical  Pulse 
Compressor.  The  power  output  from  the  3695  is  approximately  1  W.  The  pulse  is  then  sent  through  the 
Model  3500  rhodamine  dye  laser  which  is  synchronously  pumped  by  the  3800.  The  output  from  the  3500 
is  from  3  ps  to  300  fs  at  570  nm  and  150  mW.  This  pulse  duration  is  verified  with  a  Model  409  fast  scan 
autocorrelator.  The  pulses  from  the  3500  are  chirped  in  a  second  fiber  in  the  CPA-2  Chirper.  The  output 
from  the  chirper  is  in  the  nJ  range  at  570  nm  with  a  30  nm  bandwidth.  The  chirped  pulse  is  then  delayed 
and  amplified  approximately  0.5  million  times  to  100  pJ  in  a  PDA-1  three  stage  Kiton  red  dye  amplifier. 
The  dye  amplifier  is  pumped  by  the  GCR-3RA  Nd:YAG  Regenerative  Amplifier  that  operates  at  either  5 
ns  or  80  ps,  532  nm,  50  mJ,  and  an  adjustable  pulse  repetition  frequency  of  single  pulse  mode  to  10  Hz. 
The  GCR  is  seetted  by  the  3800  to  maintain  timing  of  the  pulses.  The  optical  delay  is  included  to  allow 
both  the  pulse  from  the  GCR  and  the  chirped  pulse  from  the  3800  to  enter  the  dye  amplifier  at  the  same 
time.  The  output  pulse  is  then  rephased  in  the  CPA-2  Prism  Rephaser  Module.  The  optimum  output  of 
the  laser  system  is  a  570  nm,  90  fs,  100  pJ  pulse  that  operates  in  single  shot  to  10  Hz  mode.  The  pulse 
duration  of  the  picosecond  and  femtosecond  pulses  are  measured  at  the  output  with  an  Inrad  Slow  Scan 
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Autocorrelator.  The  pulse  profile  is  approximately  80%  Gaussian  and  is  monitored  with  a  Cohu  CCD 
camera  input  to  the  Photon  BeamGrabber  beam  profiler  continuously  during  the  experiment. 

The  most  useful  feature  of  the  laser  system  is  the  wide  variety  of  wavelengths,  pulse  durations 
and  energy  levels  available  with  relative  ease  of  reconfiguration.  This  is  done  by  picking  the  pulse  off  at 
different  points  in  its  progression  through  the  system.  For  example,  300  fs  pulses  may  be  created  by 
bypassing  the  CPA-2  pulse  compressor.  In  addition,  the  laser  system  is  modular,  and  only  the  shortest 
pulses  require  all  of  the  units  to  operate  simultaneously.  All  output  pulses  have  energies  greater  than  1 00 
pj,  and  at  some  wavelengths  pulses  are  of  the  order  of  millijoules.  The  output  wavelength  is  determined 
by  the  pulse  dye  amplifier  bandwidth,  and  the  chirped  pulse  fi’om  the  dye  laser  must  match  this  bandwidth 
in  order  to  be  completely  amplified.  Otherwise,  when  the  amplified  pulses  are  recombined,  the  pulsewidth 
will  not  be  90  fs  or  shorter.  Pulses  at  a  wavelength  of  580  nm  are  obtained  for  pulsewidths  of  90,  300,  and 
800  fs,  and  1  to  3  ps.  Pulses  at  a  wavelength  of  532  nm  are  obtained  for  pulsewidths  of  50  ps  with  the 
seeded  Nd;YAG  regenerative  amplifier  and  5  ns  with  the  unamplified  Nd:YAG  Q-switched  laser,  with 
pulse  energies  for  both  in  the  millijoule  range.  If  average  power  is  needed  rather  than  energy,  modelocked 
output  pulses  at  82  MHz  are  obtained  at  532  nm  and  800  mW,  and  at  570-640  nm  at  200  mW  without 
amplification.*^ 


Figure  1:  Schematic  of  Ultrashort  Pulse  Laser  System 
used  for  LIB  threshold  measurements  at  Brooks  AFB. 


The  nanosecond  pulses  were  produced  with  a  Spectra-Riysics  DCR-11  Nd;YAG  Laser  that 
operates  at  1064  nm,  10  ns,  50  mJ,  and  from  single  shot  to  10  Hz  mode.  The  pulse  duration  is  verified 
with  a  Scientech  PIN  diode,  whose  response  curve  is  read  off  a  Tektronix  602A  digitizing  oscilloscope 
continuously  during  the  experiment.  The  pulse  profile  is  approximately  80%  Gaussian  as  verified  with  a 
Cohu  CCD  camera  input  to  Beamcode  beam  analysis  software  by  Big  Sky  Software  Corporation. 
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4.2.  Optical  Setup 


The  optical  setup  is  shown  in  Figure  2,  The  pulse  energy  from  either  the  Ultrashort  Pulse  Laser 
System  or  the  Spectra-Physics  OCR- 11  is  measured  at  A1  to  verify  the  energy  entering  the  optical  setup 
and  the  percentage  of  light  split  by  the  first  beamsplitter,  BSl.  The  amount  of  energy  split  by  a 
beamsplitter  varies  slightly  with  polarization  and  wavelength  of  the  pulse.  The  detectors  were  either  J50, 
J25,  or  J4  Molectron  detectors  output  to  a  Molectron  JD2000  Joulemeter  Ratiometer.  M2  is  placed  in  the 
setup  to  measure  the  energy  ratio  B1/A2.  During  the  experiment  the  pulse  proceeds  through  LI ,  a  17  mm 
aspheric  lens  which  focuses  die  pulse  down  to  a  spot  size  diameter  of  approximately  25  fim  in  a  cuvette  of 
water.  The  17  mm  aspheric  models  the  17.1  reduced  focus  of  the  eye*^.  The  spot  size  was  verified  with  a 
knife  edge  measurement  in  both  air  and  a  cuvette  of  water.  The  cuvette  is  an  L-shaped  jacketed  cell  that 
has  two  chambers,  one  for  the  media  and  one  connected  by  Tygon  tubing  to  the  pump  of  the  MGW  Lauda 
RC-6  temperature  controlling  bath.  The  output  pulse  from  the  cuvette  is  recollimated  by  either  another  17 
mm  aspheric  lens  or  a  25  mm  spherical  lens  with  wider  aperture  than  the  17  mm  lens.  The  energy  is 
measured  out  of  the  cell,  A3,  at  the  beginning  of  each  run  to  measure  the  ratio  of  energy  split  by  BS2. 
The  output  energy  is  recorded  at  A4.  The  other  arm  of  the  split  pulse  is  sent  into  the  Cohu  camera  for 
beam  profiling.  The  neutral  density  filters  (ND)  block  the  light  to  a  level  that  prevents  the  camera  from 
saturating.  The  plasma  is  imaged  with  two  lenses  of  focal  length  230  mm  and  1 40  mm  onto  the  Hitachi 
CCTV  camera.  The  filter  FI  is  one  of  two  beamsplitters  which  splits  either  580  nm  or  532  nm  at  45°  and 
passes  all  other  wavelengths.  The  image  is  monitored  on  a  black  and  white  TV  which  is  viewed  to 
determine  if  a  plasma  has  formed. 

The  determination  of  whether  LIB  occurred,  for  all  pulse  durations,  was  visual.  In  the 
nanosecond  time  regime,  the  visible  flash  from  the  infrared  laser  pulse,  as  well  as  the  loud  snapping  sound 
from  the  resultant  supersonic  shock  wave  made  LIB  discernment  obvious.  In  the  pico-  and  femtosecond 
time  regimes,  the  image  viewed  on  video  of  light  passing  through  the  filter  FI ,  was  sufficient  to  determine 
the  existence  of  the  LEB  event  The  gain  was  high  enough  to  saturate  the  camera  when  LIB  occurred. 
However,  in  the  pico-  and  femtosecond  regimes,  there  was  no  auditory  signal  accompanying  the  LIB  to 
aid  in  distinguishment  In  addition,  the  small  size  of  the  plasma  in  the  femtosecond  regime  made 
determination  particularly  difficult. 

Imaging  TV  BeamGiabber 


Figure  2:  Optical  setup  used  for  temperature  dependence  studies. 
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4.3.  Energy  Measurement 


The  energy  is  measured  and  recorded  at  a  number  of  points  in  the  setup  at  the  beginning  of  each 
run  to  determine  the  losses  from  particular  optical  elements.  The  energy  from  the  Ultrashort  Pulse  Laser 
System  or  the  Specnra-Physics  Nd:YAG  laser  into  the  setup  is  measured  at  Al.  The  energy  ratios  A3/B1 
and  A4/B1  are  measured  to  verify  the  losses  from  the  lenses  and  second  beamsplitter  BS2.  The  ratio  of 
energy  that  enters  into  the  lens/cell  configuration  is  measured  by  A2/B1. 

During  each  temperature  run,  the  energies  at  A4  and  B 1  are  measured  by  the  JD2000  and  sent  to 
the  National  Instruments  Labview  software  virtual  instrument  (VI).  The  ratio  A2/B1  is  multiplied  by  the 
measured  energy  from  B 1  to  get  the  adjusted  energy  into  the  cell  (E^^).  A  probability  versus  Ej„  curve  was 
generated  with  this  procedure. 


4.4.  Temperature  Control 

Temperature  control  is  accomplished  in  the  setup  with  a  jacketed  cuvette  and  temperature 
controlled  bath.  The  temperature  control  for  the  bath  is  within  ±  0.1°C.  The  temperature  control  of  the 
media  in  the  cuvette  is  within  ±  0.25°C.  The  temperature  of  the  media  is  measured  with  a  Maikson 
Model  95  Digital  Temperature  Meter  and  a  Yellow  Springs  Instrument  Model  700  probe.  The  range  of 
temperatures  measured  is  fro.m  20°C  to  80°C  in  increments  of  10^.  To  maintain  physiological 
characteristics  we  chose  not  to  measure  sample  temperatures  near  the  freezing  or  boiling  points. 


4.5.  Statistical  Analysis 

Due  to  the  probabilistic  nature  of  LIB,  a  method  of  calculating  the  thresholds  for  breakdown 
based  upon  a  normal  distribution  of  data  points  is  needed.  Probit  analysis'*  is  one  means  of  calculating 
the  probability  of  breakdown  versus  energy  dose  delivered  to  a  cuvette  of  media.  Probit  curves  and 
fiducial  limits  were  computed  on  Statistical  Analysis  System  (SAS).  Probit  curves  give  the  probability 
that  an  energy  dose  will  cause  breakdown.  The  probability  is  usually  reported  as  an  ED^^^  value  which 
gives  the  xx%  probability  that  a  certain  event  will  occur.  In  addition,  probit  analysis  gives  fiducial  limits 
of  the  probability  which  express  the  confidence  limits  of  a  probability.  Figure  3  gives  a  typical  probit 
curve  with  95%  fiducial  limits.  For  example.  Figure  3  shows  an  EDj^  of  0.275  mJ  with  upper  and  lower 
fiducial  limits  of  0.312  mJ  and  0.243  mJ  for  tap  water  at  50°C  and  laser  pulse  duration  of  60  ps.  The 
negative  energy  values  are  an  artifact  of  Probit  Analysis. 

The  error  bars  of  the  EDjq  versus  temperature  curve  are  calculated  by  taking  the  root-sum-square 
of  the  95%  fiducial  limits  and  the  instrument  uncertainty  of  the  adjusted  input  energy.  The  error  bars 
presented  in  Figures  4  and  5  may  not  fully  represent  the  uncertainty  of  the  data  for  this  size  data  set.  For 
example,  the  error  bars  should  include  any  fluctuations  in  the  laser  system  but  do  not  include  fluctuations 
in  spot  size,  which  are  impossible  to  measure  during  the  experiment.  Therefore,  the  actual  uncertainties 
may  be  somewhat  larger  than  the  ones  presented  in  Figures  4  and  5. 
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Probability 


Figure  3:  Shown  is  a  typical  probit  curve  (solid  line)  and  95%  flducial  limits  (dotted  lines) 
for  tap  water  at  50^  and  laser  pulse  duration  of  60  ps. 

5.  RESULTS 

5.1  First  Method 

The  first  method  used  to  attempt  to  resolve  a  temperature  dependence  was  described  above.  The 
temperature  was  held  constant  and  the  energy  was  varied  from  the  energy  at  which  there  is  0%  probability 
of  breakdown  to  the  energy  at  which  there  is  99%  probability  of  breakdown.  The  order  of  the  temperature 
data  was  random,  insuring  no  time  dependence  or  other  unmeasurable  effect  of  the  laser  system  that  may 
skew  the  data  over  the  course  of  the  experiment.  Approximately  200  single  shot  pulses  were  delivered  to 
the  cuvette  to  produce  a  probit  curve. 


5.1.1  Nanosecond  Time  Regime 

Figure  4  shows  the  nanosecond  EDjg  breakdown  threshold  versus  temperature.  The  difference  in 
threshold  levels  for  different  media  is  noted  and  has  been  studied  previously’^-’®  '^.  The  average  threshold 
for  bovine  vitreous  decreased  constantly  from  3.85  mJ  at  20°C  to  2.87  mJ  at  80°C.  However,  all  the  values 
fall  within  the  error  bars.  The  threshold  for  saline,  high  purity  water,  and  tap  water  exhibited  large 
variation  but  no  clear  trend.  The  large  discontinuity  in  threshold  for  high  purity  water  from  3(yC  to  40°C 
suggests  a  statistically  significant  downward  temperature  trend.  The  slightly  smaller  drop  from  40°C  to 
50^  and  the  increase  from  50°C  to  60°C  for  tap  water  also  suggest  a  statistically  significant  change  due 
to  temperature.  However,  these  observations  are  inconsistent  with  the  vitreous  and  saline  data. 
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Figure  4:  Graph  represents  the  EDj,  LIB  thresholds  vs.  temperature  for  10  ns 
laser  pulses.  The  Error  Bars  are  the  root-sum-square  of  95%  iiducial 
limits  and  uncertainty  of  energy  measurement 


5.1.2  Picosecond  Time  Regime 


Figure  5  shows  the  picosecond  EDjq  breakdown  threshold  versus  temperature.  The  results  are 
similar  to  those  found  for  the  nanosecond  time  regime.  For  saline,  all  data  points  fell  within  the 
uncertainty.  Tap  1  and  tap  2  are  runs  from  two  successive  days  and  exhibit  the  large  variation  found  in 
this  study.  The  difference  between  threshold  levels  for  equivalent  samples  may  be  due  to  different 
impurity  species  that  reside  in  each  sample.  With  the  exception  of  the  data  point  at  80°C,  all  threshold 
values  for  tap  1  fell  within  the  uncertainty.  For  tap  2,  all  data  points  except  30°C  and  IQPC  fell  within  the 
uncertainty.  As  in  the  nanosecond  data,  these  observations  are  inconsistent  and  further  smdy  is  needed  to 
explain  this  anomalous  behavior. 


44  /SPIE  Vol.  2134A  Laser-Tissue  Interaction  V  (1994) 


10  20  30  40  50  60  70  80  90 


Temperature  (C) 

Figure  4:  Graph  represents  the  EDjg  LIB  thresholds  vs.  temperature  tor  60  ps  laser 
pulses.  The  Error  Bars  are  the  root-sum-square  of  95%  fiduci^il  limits  and 
uncertainty  of  energy  measurement. 

5.1.3  Femtosecond  Time  Regime 


The  femtosecond  time  regime  yielded  the  least  quantitative  results  but  the  most  interesting 
observations  of  the  nature  of  laser  induced  breakdown.  Previous  studies  have  found  the  EDjq  threshold  erf 
LIB  to  be  approximately  1 .5  pJ  at  580  nm  and  90  fs’.  Our  study  observed  the  threshold  for  LIB  at  90  fs  to 
be  much  greater  than  1 .5  pj.  A  more  quantitative  threshold  for  LIB  at  this  pulse  duration  could  not  be 
found  because  the  maximum  energy  the  laser  can  output  is  approximately  100  pJ.  The  maximum  energy 
the  laser  can  delivered  to  the  cuvette  is  approximately  50  pJ  accounting  for  the  optics  in  the  system  and 
setup.  Therefore,  a  full  probit  curve  could  not  be  produced  due  to  the  limited  energy  output.  Another 
observation  exclusive  to  the  femtosecond  regime  is  the  decreased  recovery  time  of  the  media.  With  the 
pulses  delivered  single  shot  at  a  constant  energy  level,  no  LIB  events  were  observed  to  occur.  However,  at 
the  same  constant  energy,  many  LIB  events  were  observed  with  the  pulse  repetition  frequency  of  the  laser 
at  10  Hz.  This  decreased  recovery  time  was  not  due  to  electron  recombination,  acoustic,  or  thermal 
effects,  which  occur  on  the  order  of  approximately  lO'^  10^,  and  lO"^  seconds  after  the  LIB  event, 
respectively.  Gas  and  bubble  formation  may  explain  the  observed  response  of  the  media  to  femtosecond 
pulses  at  10  Hz. 

In  addition,  self-phase  modulation  (SPM)  and  bubble  formation  were  observed  to  occur  at 
approximately  1-2  pJ.  This  result  is  contrary  to  the  results  found  in  previous  studies  for  longer  pulse 
durations*®’*®  where  the  threshold  for  SPM  is  very  close  to  the  threshold  for  LIB.  This  lends  further 
credence  to  the  postulate  of  two  distinct  types  of  plasmas  at  this  pulse  duration  each  at  a  different 
irradiance  threshold**.  The  hot  plasma,  occurring  at  the  higher  threshold,  produces  a  visible  broadband 
flash  in  all  directions  and  is  similar  to  the  plasmas  produced  at  other  durations  with  the  exception  of  the 
threshold  level.  The  cold  plasma,  occurring  at  a  much  lower  threshold,  is  without  a  broadband  flash  but 
with  other  characteristics  associated  with  LIB,  namely  initiation  by  multiphoton  ionization,  continuum 
generation,  and  cavitation  bubble  formation.  Continued  study  is  needed  to  fully  understand  LIB  at  these 
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ultrashort  pulse  durations. 


5J1  Second  Method 

The  inconsistencies  presented  above  in  the  first  method  led  us  to  seek  more  precise  experimental 
techniques  to  measure  any  temperature  dependence.  A  second  method  was  devised,  in  which  the  energy 
was  held  constant  somewhere  within  the  0%  and  99%  probability  of  breakdown  while  the  temperature  was 
varied.  Approximately  200  single  shot  pulses  were  delivered  to  the  cuvette  of  media  at  a  constant  energy 
and  temperature.  The  energy  was  windowed  to  a  set  energy  range  and  for  all  cases  had  a  standard 
deviation  less  than  one  tenth  of  the  average  energy  delivered.  Probit  analysis  was  not  used,  with  the 
expectation  that  by  reducing  as  many  areas  of  probability  as  possible,  a  more  realistic  threshold  could  be 
obtained  and  thus  any  slight  variation  in  the  temperature  dependence  measured. 

Table  1  lists  the  percentage  of  negative  responses  obtained  for  uiree  temperatures  for  tap  water, 
high  purity  water,  and  saline  at  constant  pulse  energy.  A  decreasing  percentage  with  increasing 
temperature  would  indicate  a  temperature  dependence.  Saline  and  tap  water  exposed  to  60  ps  pulses  show 
the  hypothesized  decrease  in  percentage.  The  other  sets  of  data  all  exhibit  a  large  variation  and  no  clear 
trend. 


Table  I:  Percentage  of  Negative  Responses  for  constant  energy  and  varying  temperature 


60  Picosecond 

10  Nanosecond 

Saline 

Saline 

24<C 

50% 

66.C 

69% 

79% 

72% 

62% 

509C 

17% 

48% 

68% 

63% 

43% 

71% 

80°C 

34% 

39% 

52% 

;  72% 

53% 

58% 

6.  DISCUSSION 


6.1  Experimental  Discussion 

The  data  presented  above  show  no  clear  trend  of  temperature  dependence.  Some  of  the  data  may 
appear  to  exhibit  a  trend,  such  as  the  bovine  vitreous  in  the  first  method  and  the  saline  in  the  second 
method.  However,  the  large  variability  of  all  the  thresholds  discount  any  attempt  to  reach  a  positive 
conclusion  based  upon  this  data.  This  was  due  to  the  difficulty  in  consistently  measuring  the  thresholds 
which  are  based  on  a  probability  curve  without  taking  an  inordinately  large  sample  size.  There  were  a 
multitude  of  effects  that  needed  to  be  controlled  to  differentiate  any  supposed  temperature  dependence 
from  the  other  factors. 

6.2.  Theoretical  Discussion 

Our  measurements  on  ocular  media  using  various  pulsewidths  showed  no  consistent  temperature 
dependence  large  enough  to  resolve  with  current  experimental  techniques  in  the  range  from  20“C  to  80°C. 
A  first-order  theoretical  analysis,  using  optical  absorption  data  for  water,  seems  to  support  the  conclusion 
that  temperature  dependence  over  this  range  will  be  small  to  negligible. 

There  are  two  ways  increased  temperature  could,  in  theory,  reduce  the  optical  breakdown 
threshold  of  water  and  aqueous  solutions.  First,  higher  temperatures  would  increase  the  probability  of 
molecular  ionization  through  thermal  collisions.  If  the  increase  in  thermal  ionization  was  significant  over 
the  temperature  range  considered,  then  the  avalanche  breakdown  threshold  would  decrease  due  to 
increased  initial  (seed)  electron  density.  This  effect  is  unlikely  due  to  the  large  ionization  energy  of  water. 
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approximately  6.5  eV“'^'.  At  80°C,  kgT  -  0.03  eV,  a  thermal  energy  far  too  low  to  give  a  statistically 
significant  probability  of  ionization. 


The  second  possibility  is  that  higher  temperatures  could  reduce  the  average  ionization  energy  of 
the  medium,  which  would  reduce  the  threshold  for  both  avalanche  and  multiphoton  breakdown.  UV 
absorption  measurements  in  water,  however,  show  that  for  temperature  variation  from  24°C  to  80°C, 
transition  energy  decreases  only  0.12  eV,  roughly  a  2%  decrease^.  Such  a  small  decrease  is  unlikely  to 
produce  measurable  changes  in  the  breakdown  threshold. 

The  comments  above  exclusively  apply  to  ionization  of  water  molecules.  It  is  possible  that 
impurities  with  very  low  binding  energies  may  create  a  temperature  dependent  breakdown  in  aqueous 
solutions.  Measurements  in  impure  media  such  as  tap  water,  however,  have  not  shown  this  effect. 

6.3.  Conclusion 

The  conclusion  of  this  study  is  that  the  LIB  measured  thresholds  are  not  dependent  on  the 
temperature  of  the  media  in  the  range  studied.  Furthermore,  the  LIB  threshold  data  measured  at  21®C  to 
25°C  presented  in  literature  are  valid  comparisons  to  the  LIB  thresholds  that  occur  at  35°C  to  37°C  in  the 
eye. 

In  addition,  valuable  collateral  information  has  been  gained  from  this  study.  The  observations  of 
femtosecond  pulses  presented  above  in  particular  have  led  to  discussion  of  the  physical  properties  of  the 
plasmas.  The  differences  in  the  breakdown  threshold  and  superbroadening  threshold  at  this  pulse 
duration  could  indicate  an  occurrence  of  nonlinear  processes  and  need  to  be  studied  further.  The 
decreased  recovery  time  of  the  media  in  the  femtosecond  time  regime  is  also  puzzling  considering  most 
side  effects  of  LIB  have  ceased  by  milliseconds  after  the  event.  In  addition,  the  femtosecond  breakdown 
thresholds  need  to  be  reevaluated  to  either  validify  or  reftite  earlier  threshold  measurements. 
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